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Abstract We investigate the interacting agegraphic dark 
energy in Brans-Dicke theory and introduce a new se- 
ries general forms of dark sector coupling. As examples, 
we select three cases involving a linear interaction form 
(Model I) and two nonlinear interaction form (Model 
II and Model III). Our conclusions show that the ac- 
celerated scaling attractor solutions do exist in these 
models. We also find that these interacting agegraphic 
dark energy modes are consistent with the observational 
data. The difference in these models is that nonlinear 
interaction forms give more approached evolution to the 
standard vlCDM model than the linear one. Our work 
implies that the nonlinear interaction forms should be 
payed more attention. 

PACS 95.36. +x ■ 98.80.-k • 98.80.Es 
1 Introduction 

Recent astronomical observations have provided strong 
evidence that our universe is undergoing an accelerated 
expansion due to an exotic energy component with neg- 
ative pressure which is called dark energy [UEEIl] ■ By 
far a leading candidate of dark energy model is the Ein- 
stein's cosmological constant model (ACDM) which is 
consistent with the late-time observational data. How- 
ever, it is confronted with the so-called "cosmologi- 
cal constant problem" and "coincidence problem" [5]. 
Therefore, many dynamical dark energy models [5], such 
as quintessence [7], phantom [S], quintom [S], tachyon 
[TO] , generalized Chaplygin gas [TTj, etc, have also been 
taken into account. 

The cosmological constant problem may be essen- 
tially an issue of quantum gravity problem, since the 
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cosmological constant is commonly considered as the 
vacuum expectation value of some quantum fields. Al- 
though a complete theory of quantum gravity has not 
been established, according to some principles of quan- 
tum gravity, one can make some attempts to probe the 
nature of dark energy model. The holographic dark en- 
ergy model is just an appropriate and interesting exam- 
ple. This model is based on the holographic principle of 
quantum gravity theory [T2], and is derived from the 
relationship between the ultraviolet (UV) and the in- 
frared (IR) cutoffs proposed by Cohen et al. in ref. [13] . 
According to the limit set of the formation of a black 
hole, the UV-IR relationship gives an upper bound on 
the zero point energy density p q = 3L~ 2 /8irG, which 
means that the maximum entropy of the system is of 
the order of S^g , where L is the scale of IR cutoff and 
Sbh is the entropy of the black hole. Compared with 
the holographic dark energy model and based on the 
Kdrolyhdzy relation [Tl] , the so-called agegraphic dark 
energy model was proposed, where the age of the uni- 
verse T = J dt is used as the IR cutoff L [15] . Further- 
more, in ref. |16j , the interacting agegraphic dark energy 
has been introduced and investigated. It was shown that 
the equation of state of interacting agegraphic dark en- 
ergy can cross the phantom division. The interacting 
agegraphic dark energy model also has been extended 
to the universe with spatial curvature in ref. [17] . Re- 
cently, using the phase space analysis, it was shown that 
the accelerated scaling attractor solutions of the inter- 
acting agegraphic dark energy in the Einstein universe 
did exist and the results agree with the observations 
PI- 

On the other side, scalar-tensor theories of gravity 
have been widely applied to cosmology [T5]. The sim- 
plest alternative to Einstein's general relativity which 
includes a scalar field in addition to the tensor field 
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is Brans-Dickc theory. This theory is more consistent 
with the Mach's principle and less reliant on absolute 
properties of space |20j . It got a new impetus in re- 
cent years because it arises naturally as the low energy 
limit of many theories of quantum gravity such as super 
string theory or Kaluza-Klein theory. Noticing that the 
holographic dark energy density belongs to a dynamical 
cosmology constant, it is more natural for a dynamical 
frame to replace the Einstein's general gravity. There- 
fore, it is worthwhile to investigate the holographic dark 
energy model within the framework of Brans-Dicke the- 
ory [2T1I221I231I24] . The extended holographic dark en- 
ergy model with Hubble horizon in Brans-Dicke theory 
has been proposed and it is found that the model is 
not a viable dark energy model unless the Brans-Dickc 
scalar field has a potential. So it is a very interesting at- 
tempt to deeply investigate the agegraphic dark energy 
model in the framework of Brans-Dicke theory . 

Considering that dark energy (DE) and dark matter 
(DM) contribute to the most fraction of the content of 
the universe, it is natural to look into the possibility of 
the interaction between DE and DM, which has been 
widely discussed It has been argued that the cou- 
pling between DE and DM can provide a mechanism 
to alleviate the coincidence problem and lead to an ac- 
celerated scaling attractor solution with similar energy 
densities in the dark sector today [2"M2"Tll28[|29j . Notic- 
ing that there is no fundamental theory which can be 
used to select a specific interacting dark energy model, 
any interacting dark energy model will necessarily be 
phenomenological. There are two criterions to deter- 
mine whether the model is correct and feasible. One is 
the observations, the other is to examine whether the 
interacting model can lead to the accelerated scaling 
attractor solutions, which is a decisive way to achieve 
similar energy densities in dark sector and alleviate the 
coincidence problem. In this work, we firstly introduce a 
new series of interacting agegraphic dark energy models 
including linear and nonlinear forms. Using the phase- 
plane analysis, it is found that the accelerated scaling 
attractor solutions do exist in these models. What's 
more, these agegraphic dark energy models are in ac- 
cordance with the late-time observational data. 

Our paper is organized as follows. In section [21 the 
agegraphic dark energy models in Brans-Dicke theory is 
constructed and a series of dark sector coupling forms 
are introduced. In section[3J using the phase-plane anal- 
ysis, the accelerated scaling attractor solutions are dis- 
cussed in these models. In section^ using the newly re- 
leased Hubble parameter data [3Tflnn"ll3"2"ll3"3"] , these age- 
graphic dark energy models are tested. Some conclu- 
sions will be presented in section [5j 



2 Construction of the interacting agegraphic 
dark energy in Brans-Dicke theory 

According to our metric convention, (+,—,—,—), the 
Lagrangian for Brans-Dicke theory with a scalar field 
in the Jordan frame is 

1 



Cbd = y/=g[(-<pR + u-grdn<pd v v) +£ M (#)], (1) 
V 

where the dimensionless u> is the coupling constant, and 
is the matter Lagrangian. The current obser- 
vational constraint on u is u> > 10 4 ; it recovers Ein- 
stein's general relativity when u — > oo. In particular, it 
is expected that (p(t, x) is spatially uniform and evolves 
slowly only with cosmic time t so that (p(t,x) — > <p(t). 
Following ref. |35j . we can introduce a new field <j) as 
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So the Lagrangian for Brans-Dickc theory can be writ- 
ten as following form 
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,<j> + C M (p)\. (3) 



Considering a classical perfect fluid with the energy- 
momentum tensor = diag(p, — p, — p, — p), the grav- 
itational field equations derived from the variation of 
the action Eq. ^ with respect to the flat Robertson- 
Walker metric are 

3 9 " 9 1 '.9 . 3 • 

-H<jxp = p, 



4u> 
1 

4u/ 



l H 2 



2 2uj 
2 {2H + 3H 2 
3 

P) 



1 



##-( 



1 



2uj' 



1 

— c 

2uj 



3H<j> 
3H(p 



-(H- 
0. 



2H 2 



where the dot is the derivative with respect to time and 
H = ^ is the Hubble parameter. Combining the above 
equations, we get 
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[<jjp+(w+2)p]. 
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The total matter is supposed to be composed of two 
parts: pressurclcss dark matter and agegraphic dark en- 
ergy. So the total energy density p includes the energy 
density of agegraphic dark energy p q and the energy 
density of dark matter p m in this model. In Brans-Dickc 
theory, the effective gravitational constant G e ff can be 
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naturally be defined as 
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the agegraphic dark energy can 
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where n is a positive constant. Now we turn to con- 
sidering the interaction between dark matter and dark 
energy. The balance equations of the agegraphic dark 
energy and dark matter can be written respectively as 



p q + 3H(l + ui q )p q 
Pm + 3Hp m = Q, 



-Q, 



(10) 

(ii) 



where Q denotes the phenomenological interaction term. 
The interacting term is always a function of the Hubble 
parameter H, the density of dark energy, and the den- 
sity of the dark matter. One can often find the following 
interaction model forms [I51|T5pSP5P7P5pgpnp51 



Q = 3Hap q , 3H(3p m , 3jH(p q + p m ), 



3H(ap q + /3p m ), 3/3Hp°!p, 



q r rn J 



(12) 



where a, /?, 7 are positive constant parameters. We would 
like to consider the more general interaction as 
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Q = 3Hp q g{t) = 3Hp q ]T M\ 



(13) 



where £ = ^ , g(£) = A i^-> and A i is a positive 

constant parameter with respect to £\ Then Eq. lfTU)) 
and Ea. ([TT|) can be written as 

p q + 3H{l + u q )p q = -3H Pq g(0, (14) 
p m + 3Hp m = 3Hp q g{£). (15) 

In this paper, as examples, <?(£) will be selected as: 

Model I: gi(0=a 1 +M, (16) 
Model II: ff2 (0 = a 2 +hi 2 , (17) 
Model III: g 3 (0 = a^" 1 + M- (18) 

It is obvious that Model I is linear interaction, while 
Model II and III are nonlinear interactions, which can 
be taken as the extending of Eq. (fl"2"]) . 

3 Phase-space analysis 

According to the discussion in section[2j it is easy to see 
that cqs. fll]), ©, flSJ), (H]), and (TJS]) can give a closed 
system which can determine the cosmic behavior. In 
order to study the dynamical behavior of interacting 
dark energy and dark matter, we further introduce the 
following dimcnsionless variables 
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wip pm > x= w (19) 



Similar forms of selecting these dimcnsionless variables 
in agegraphic dark energy model was firstly used in ref. 



[T5] . Using these dimensionlcss variables, Eq. be- 
comes 
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Now, the critical density can be defined as p c — g7rG - — 

'^f- gl], so Eq. ([201) becomes fl q + Q m + [2$ = 1, 
where 



f2 q — x , firn y , f2^ 
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From Eq. (f^Tj) . we find that the Brans-Dicke scalar field 
plays a role of dark energy, so we can assume that 
both the agegraphic and the scalar field drive our uni- 
verse to accelerate [24] . So we have 
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Noticing that < f2 q , f2^, fi m < 1, the limit of x, y, 
and A can be obtained as0<a;<l,0<y<l, 

*.( 1 - y fi+% ) < A < Oor 1 < A < £(1- yiTf ) 
respectively. 

Subsequently, thinking of Eqs. flSJ), we have 

x 



f(x,y,\) 
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A' = — 3A - A 2 H h ( A) f(x, y, A), 



(24) 
(25) 
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where x' = dx/dN,y' = dy/dN,X' = dX/dN,N 

2 

lna,£ = and 



(26) 



/(,, y ,A) = ^=4A-2^-^ 

The state parameter w g for the agegraphic dark energy 
could be expressed in terms of these new variables as 



(27) 
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For completeness, we give the deceleration parameter 

-l-f(x,y,X). (29) 
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In the interacting agegraphic dark energy models, the 
properties of agegraphic dark energy are determined 
by the parameters Aj, u> and n. Eqs. ((2"4")) - ([2"T|) are just 
the functions of x,y, and A, not the functions of N or 
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other variables. So this dynamical system is just a three- 
dimensional autonomous system. For an autonomous 
system X' = /(X), there are some critical points X c 
satisfying X' = 0, so we have 

c 



x c (f{x c ,ycAc) + 



-3A C — A,, 



-) = o, 

11 

) + A c + f(x c ,y c , A c )] 



Vc 



,y c ,X c ) = 0, 



(30) 
(31) 
(32) 



where £ c = %. In order to determine the stability prop- 

x c 

erty of the critical points, it is necessary to expand 
the autonomous system X' = /(X) around the crit- 
ical points X c . Setting X = X c + U, where U is a 
column vector of the perturbation of the variables, one 
can expand the equation for the perturbation up to the 
first order as U' = M • U, where the matrix M contains 
the coefficients of the perturbation equations. The sta- 
bility property for each critical point is determined by 
the eigenvalues of M [25]. The 3x3 matrix M of the 
linearized perturbation equations is 
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We can examine the sign of the real part of the eigen- 
values of M-x=x a ,y=y a ,\=\ a j which determines the type 
and stability of the critical points (x c , y c , X c ). The sim- 
plest finite critical points and their properties for this 
model are summarized in the Table [TJ Here we only con- 
centrate on the attractor, so the real part of the eigen- 
values of 'M- X =x,,y=y,,\=x, should be negative. If the 
autonomous system of Eas.([24"|). ([25]) and ((26)) presents 
scaling solutions, the coincidence problem gets substan- 
tially alleviated because, regardless of the initial condi- 
tions, the system evolves toward a final state where the 




Fig. 1 The phase plane for the model I with n = 20, lo = 10 3 , 
ai = 0.3 and /3i = 0.16. The red star stands for the late-time 
attractor with x, « 0.86, y, w 0.26, \, « 9.9 x 10~ 4 . 



ratio of dark matter to dark enenrgy stays constant. 
Using Eq. ([2"3"]) , it is easy to find that 



-(1 



r 2 )[2xx' + (— A-2)A'] 
3 



(33) 



The scaling solutions mean r' = 0. Obviously the at- 
tractor (a;*, j/*, A*) is just a scaling solution, because 
x'\x, =0,A'|a„ = lead to r'\ Xry y r ^\ r — 0. What's more, 
if the critical point is the accelerated scaling attrac- 
tor solution, the following constraints q(x*, y*, A*) < 
0, r(x*, y», A*) ss ^|| « 0.353 should also be satisfied. 

Obviously the accelerated scaling attractor solutions 
can be obtained by fixing the parameters Ai, ui, and n. 
We have found the accelerated scaling attractor solu- 
tions in Model I, Model II, and Model III using the nu- 
merical analysis method. The results are shown in Fig. 
[U f2J [3] respectively. The numerical solutions show that 
the cosmic evolution is insensitive to the initial condi- 
tions. The different interacting model forms just affect 
the intermediate evolution. This makes the coincidence 
problem to be substantially alleviated. 

Figures [H [5j [6] show that the three interacting age- 
graphic dark energy models give very similar conclu- 
sions. From Fig. 0] one can see the cosmic evolution in 
these models, where £1$ dominates the universe at the 
early epoch and then the matter is becoming increas- 
ingly dominant over J?^ afterwards. After a transitory 
matter-dominated epoch, the universe becomes accel- 
erated expansion in the agegraphic dark energy domi- 
nated epoch. From Fig. [5] the universe is now undergo- 
ing an accelerated expansion phase. The similar curves 
of q in three models show that this cosmic acceleration 
is arising recently. The evolution of equations of state of 
dark energy is presented in Fig. [5] There is a phantom- 
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Table 1 Location of the critical points of the autonomous system of Eqs.j2l]|, Q25I) and (12611 . their stability 
and dynamical behavior of the Universe at those points. 



(^cjJ/cAc) coordinates Stability character 



(0,0, Ai = ^j(l± y/T + ^f)) unstable -4 - 4Ai 

(Q, y 2 = yjl- 1f\\ + 2A 2 , A 2 = [1 ± ^1 + 



unstable 



3-2ojX, 



(X3,0,A3 = £[1± Jl + ¥•(!- a:§)]) unstable -1 + ^ 



(x, = " (6 " A j+" A2 ~ 6) , y„ = ,/l - +2A« -a£, A,) attractor 




Fig. 2 The phase plane for the model II with n = 20, lj = 
10 3 , a 2 = 0.36 and /3 2 = 10 -6 . The red star stands for the 
late-time attractor with x* w 0.86, y„ w 0.26, A, ss 9.9 x 
10" 4 . 




Fig. 3 The phase plane for the model III with n = 20, ui = 
10 3 , a.3 = 0.12 and /?3 = 0.1. The red star stands for the late- 
time attractor with x* « 0.86, y, « 0.26, A, « 9.9 x 10 -4 . 



like to quintessence-like transition while the values of 
1. 



u> q is always around 



4 Observational data 

Now, we would like to check the three interacting age- 
graphic dark energy models using the observational data 
and compare them with the standard ylCDM model. 

The continuity eqs. (|14p and (| 15[) can also be written 
in the standard forms 



Pq 
Pm 



3H(l+u)')p q = 0, 



3H(1 + 
where u>' = — 1 



= 0, 
" I). < 



(34) 
(35) 



_ MSI 

"q ~" * 3 V n<" *~m f ■ 

The Brans-Dicke field can be assumed as a power 
law of the scale factor 1231 



= <Poa 



(36) 



Noticing that the scalar field £2$ in figure H] is very in- 
finitesimal in the late-time universe, we can assume that 
A is a constant A w A c , and the product \ c lo results or- 
der unity, which is consistent with refs. |23ll4Tj . 

Substituting Eqs. (J3J]), ([35]), and ([3BJ into Eq. g]), 
one can get the evolution of Hubble parameter H = 
HqE(z), where the expansion rate E(z) can be expressed 
as 



E\z) = (1 
[f2 m0 exp(3 



2X r 



2w , , , 

3A! } + 



2A,. 



^dz') + ^ exp(3 ri±4dz')] 
z Jo 1 ' z 



(37) 



where the subscript "0" denotes the present value. 

We will use the newly released Hubble parameter 
data |50U3TTI53U551l54l which is directly related to the 
expansion history of the universe by its definition: H = 
a I a to determine a preferable interacting model. In the 
meantime, we choose the standard ylCDM model as the 
fiducial model which fits the observations best. The re- 
sults are presented in Fig. [7J 
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Fig. 4 Behavior of Q m (black), Q q (re&) and (green) as a Fig. 5 Behavior of deceleration parameter q as a function of 
function of N = In a for the particular models N = In a for the particular models 
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Hubble data 

Model I 

— A CDM 
■ Model II 
Model III 



1.5 



Fig. 7 The theoretical evolution of H(z) with observational 
data points (black): ylCDM (red solid), Model I (blue solid), 
Model II (blue dotted), and Model III (blue dashed) 



Apparently, the nonlinear interaction forms (Model 
II and Model III) give better evolutions to the stan- 
dard ylCDM one. The linear one shows evident devia- 
tion from the /1CDM and observation data. The nonlin- 
ear interaction between dark energy and dark matter is 
better to describe the real physical process. This com- 
parison is incomplete and it is only a rough overview. It 
is necessary to further explore the nonlinear interaction 
between the dark section of the universe. 



5 Discussions and conclusions 



-0.96 - 



Model III 



Fig. 6 Behavior of effective equation of state of dark energy 
uj' q as a function of N = In a for the particular models 



Holographic dark energy model is an interesting at- 
tempt to investigate the nature of dark energy in the 
framework of quantum gravity. Considering that the 
simplest alternative to Einstein's general relativity which 
includes a scalar field in addition to the tensor field is 
Brans-Dickc theory and the holographic bound can be 
satisfied for both the k = and k = — 1 universe in 
the Brans-Dickc cosmology [32; it is natural to extend 
the research to the holographic dark energy models of 
this theory. In this paper, we have investigated the in- 
teracting agegraphic dark energy in the flat ( k = ) 
Brans-Dicke cosmology. Firstly, a series of new general 
forms of dark sector coupling are introduced and the 
accelerated scaling attractor solutions have been found. 
Moreover, using the newly released Hubble parameter 
data, we have also tested these interacting agegraphic 
dark energy models. 

The interacting term can be selected as a function 
of H, the density of dark energy, and the density of the 
dark matter. According to this requirement, we have 
proposed the general interacting agegraphic dark en- 
ergy models. Three cases including a linear interaction 
form (Model I) and two nonlinear interaction forms 
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(Model II and Model III) have been investigated. Us- 
ing the phase-plane analysis, the dynamical behavior 
of these models has been investigated and it was found 
that the accelerated scaling attractor solutions did ex- 
ist in these models. This can alleviate the coincidence 
problem. 

Afterwards, using the newly released Hubble param- 
eter data, we tested these dynamical dark energy mod- 
els. These interacting agegraphic dark energy modes 
have given a series of reasonable pictures of the cos- 
mic evolution and they were consistent with the late- 
time observational data. In particular, we found the 
nonlinear interaction forms (Model II and Model III) 
gave more approached evolution to the standard ACDM 
model than the linear one (Model I). Our work show 
that we should pay more attention to the nonlinear in- 
teraction forms rather than the linear form. This de- 
serves further investigations. 
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